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ABSTRACT f- <>V 


A method for analyzing biaxial- and shear-loaded anisotropic rectangular panels with 
centrally located circular and elliptical cutouts is presented in the present paper. The method 
is based on Lekhnitskii’s complex variable equations of plane elastostatics combined with a 
boundary collocation method and a Laurent series approximation. Results are presented for 
anisotropic panels with elliptical cutouts and subjected to combined shear and compression 
loading. The effects on the stress field of panel aspect ratio, anisotropy, cutout size, and 
cutout orientation are addressed. Angle-ply laminates, unidirectional off-axis laminates, and 
[(±45/0/90) 3 ] s , [(±45/0 2 )3]s- and [(±45/90 2 )3]s laminates are examined. 


INTRODUCTION 


Stress distributions in laminated composite panels with cutouts are an important 
consideration in aircraft design and analysis. Cutouts are often necessary in aircraft 
structures to form access ports for electrical and mechanical systems. In addition, significant 
weight savings can be achieved through the introduction of cutouts in wing ribs and other 
aircraft components. The effects of cutout shape and orientation on the magnitude and 
distribution of the stress field are important in the design of these components. Presently, the 
majority of stress analysis methods are based on classical infinite plate theory or finite 
element analysis. Finite element analysis has been a popular approach to the stress 
analysis of finite-dimensional panels with cutouts (refs. 1 , 2). These analyses produce 
accurate results, but they are costly methods of performing parametric studies in which 
several different materials and geometries must be considered. The method described in the 
present paper is an alternative to finite element analysis for panels with cutouts. 

The method presented herein is based on the complex variable equations of plane 
elastostatics presented by Lekhnitskii (ref. 3). These equations are used in conjunction with 
a boundary collocation method and a Laurent series approximation to analyze the stress 
fields of anisotropic panels with centrally located cutouts. This method allows for biaxial, 
shear, and combined loading, and accommodates elliptical cutouts of arbitrary orientation, as 
well as circular cutouts. 

An analytical study of the stress distribution in anisotropic panels with cutouts is 
presented to illustrate the versatility of the method, and some comparisons are made with 
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experimental data. The effects of panel aspect ratio, anisotropy, cutout size, and cutout 
orientation on the stress field are examined. Angle-ply laminates, unidirectional off-axis 
laminates, as well as [(±45/0/90) 3 ]s, [(±45/0 2 ) 3 ]s, and [(±45/90 2 ) 3 ]s laminates, are considered. 
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ellipse minor axis 

constant coefficients of Laurent series 

coefficients of Laurent series constant coefficients 

circular cutout diameter and ellipse major axis 

average elastic modulus in x-direction 

average elastic modulus in y-direction 

Airy stress function 

applied forces at panel boundaries 

average shear modulus 

panel length 

arc length on interior or exterior boundary 
panel width 

x-component of boundary traction 
y-component of boundary traction 
complex variable defined as z k = x + p k y 
average strain in x-direction 

coefficient of mutual influence of the second kind which 
characterizes shearing in the xy-plane caused by a normal 
stress in the x-direction 

coefficient of mutual influence of the second kind which 

characterizes shearing in the xy-plane caused by a normal 

stress in the y-direction 

fiber orientation angle 

complex roots of the characteristic equation 

Poisson’s ratio 

inclination angle of elliptical cutout 
average normal stress in the x-direction 
average normal stress in the y-direction 
average shear stress 

functions of z k which make up the stress function 
first derivative of <|> k 
transformation variable 


A bar over a quantity denotes its complex conjugate. 
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ANALYSIS 


Theory 

The objective of this analysis is to develop a method for determining the stress 
distribution in a finite anisotropic panel with a centrally located elliptical cutout. The panel is 
loaded by in-plane forces which do not vary through the thickness, and a state of generalized 
plane stress is assumed. Thus, average material properties are employed in the present 
stress analysis, which is based on Lekhnitskii’s complex variable equations (ref. 3). A 
complete description of this analysis is provided in references 4 and 5, and a summary 
follows. 

The generalized biharmonic equation for an anisotropic material in terms of an Airy 
stress function, F, and the average material properties is 

3 4 F 3 4 F / IxdlF 

3y* - 2T lxy,x 8y3ax - |^ v xy - Gxy Jay 23 x 2 - ^xy.y Ey 3y3 X 3 + E y 3x4 - « 

This biharmonic equation can be simplified using the transformation 

z\ = x + Pt y z 2 = x + p 2 y 

where p-|, p 2 , and their complex conjugates are the roots of the characteristic equation 


M- 4 ' 2l lxy,x 




Using the above transformation, the generalized biharmonic equation can be written as 


_3_ _3_ _d = _ 
dz 2 3 z 1 


£ 

3 Z 2 


F = 0 


The solution for the stress function, F, is 


f = pi ( z i ) + <M z 2) + <M z i) + ^ 2 (^ 2 ) 


Applied displacements or tractions along the panel boundary can be related to the complex- 
valued stress function. In this paper, only traction boundary conditions are considered. The 
boundary tractions in the x-direction, X n , and in the y-direction,Y n , are related to the stress 
function by 
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where S is the length of the boundary arc that begins at a point z 0 and ends at z, and O k (z k ) is 
defined as 




If the value of this stress function is known for every point within the panel boundaries, the 
stress distribution in the panel can be determined. In the present analysis, this stress function 
is represented by a truncated Laurent series containing unknown constant coefficients: 


^k) =£B kn z£ 


N 


The resultant force on every arc of the interior and exterior boundaries is known and, 
therefore, Lekhnitskii’s force equations can be used to solve for the unknown coefficients 
through boundary collocation. Satisfying the force loading conditions along the interior and 
exterior panel boundaries results in a system of equations that can be solved for the 
unknown Laurent series constants: 


[ Cmkn ] { B kn } - { F m } 

where B kn are the unknown Laurent series constants, C mkn are the coefficients of the Lauren: 
series constants, and F m are the applied force resultants. In order to improve solution 
convergence without increasing computation time, a least squares approach is implemented. 
The least squares method allows for an increase in the number of force equations 
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considered without an increase in the number of terms in the Laurent series. In the present 
case, twice as many equations as unknowns are considered; therefore, 

[ C m kn ] is a 16N x 8N matrix 
{ Bkn } is an 8N vector 
{ F m } is a 16N vector 

To solve the system of equations for the unknown Laurent series coefficients, each side of the 
matrix equation is multiplied by the transpose of [ C m kn ]• To improve the conditioning of the 
system matrix, the following mapping function which maps all of the points inside the panel 
boundaries to the exterior of a unit circle is used: 







\ 

) 
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where D is the major axis of the elliptical cutout, and A is the minor axis. The use of this 
mapping function eliminates the problem of small numbers being raised to high powers. 
After the system of linear algebraic equations is solved for the Laurent series constants, the 
average stresses in the panel can be calculated using the following stress equations: 


crF 2 ' 2 ' 

a x = ^2 = 2Re[p 1 0 1 (z 1 ) + p 2 <D 2 (z 2 ) ] 

3 2 F 

a v = ax 2 = 2Re [ *i(zi) + ^ z 2) 1 
3 2 F 

CT *y = ' 3x3y = ' 2Re I + M> 2 ( z 2 ) ] 


Model 

A FORTRAN program was written to implement the present analysis. The 
configuration analyzed consists of a rectangular panel with a centrally located cutout as 
shown in figure 1 . The panel is of length L (in the x-direction) and of width W. The panel has 
either a circular or an elliptical cutout at its center. The circular cutout is of diameter D; the 
elliptical cutout has a major axis D and a minor axis A. The major axis of the ellipse may be 
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inclined at some angle to the y-axis or may be aligned with the y-axis. The panel is 
subjected to an applied shear stress, c xy , and/or an applied compressive stress, o x . In this 
study, when both shear and compression stresses are applied, they are equal in magnitude. 
The boundary of the cutout is assumed to be traction-free. The lamina fiber angles, 0, are 
measured with respect to the x-axis. 


RESULTS AND DISCUSSION 


Comparison with Experiment 

To assess the accuracy of the current analytical approach, comparisons are made 
between experimental data and analytical results. The stress distribution near a circular 
cutout in a rectangular orthotropic panel with a uniaxial applied compressive stress is shown 
in figure 2. The (±45/02/±45/02/±45/0/90)2s graphite-epoxy laminate has a width of 4.49 
inches, a length of 10 inches, and cutout-diameter-to-panel-width ratios, D/W, of 0.11, 0.17, 
0.22, and 0.33. Normal stress along the y-axis, a x (0,y), normalized by the applied 
compressive stress, o° is plotted as a function of the distance from the center of the panel, y, 
normalized by half the panel width, W/2. The analytical results are represented by solid lines, 
and the experimental data are represented by symbols. The experimental data were 
obtained from a study presented in reference 6. The analytical results agree well with the 
experimental data. 

A comparison between analytical and experimental data is also made for a square 
[(±45/0/90)3 ]s graphite-thermoplastic laminate subjected to an applied shear stress, a xy . The 
panel has a width of 12 inches and cutout-diameter-to-panel-width ratios, D/W, of 0.063, 
0.125, and 0.25. The normal strain along the y'-axis, e x .(0,y'), is plotted in figure 3 as a 
function of the nondimensional distance from the center of the panel, y'/W. For the purpose of 
this comparison, the panel is oriented at a 45° angle to the y-axis as shown in figure 3. The 
analytical results are represented by solid lines, and the experimental data from reference 7 
are represented by symbols. The analytical results agree reasonably well with the 
experimental data; differences between analytical and experimental results are suspected to 
be largely due to the manner in which the load is introduced. In the analysis, a uniform shear 
stress is applied to the panel boundaries. In the experiment, loads are introduced through 
the use of a picture frame test fixture, which tends to concentrate the loads at the panel 
corners. 


Analytical Study 

To demonstrate the capabilities of the analytical method, stresses are calculated for e 
range of panel dimensions, ply layups, and loading conditions. Results are presented in 
terms of a maximum stress value which is normalized by the corresponding applied stress. 
This maximum stress value is obtained by conducting a survey of stress values at numerous, 
evenly-spaced points on the panel. Points inside the cutout region are excluded. The 
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laminates analyzed are made of Hercules, Inc. AS4/3502 graphite-epoxy unidirectional tape 
and the lamina properties for the laminates studied are presented in Table 1. 

The FORTRAN program which implements the present method consists of two lines of 
input describing the panel geometry, loading conditions, and average material properties. 
Due to the minimal amount of input, the following cases require very little user time for the 
creation of the models. In addition, the program run time is short enough that the program 
can be run interactively. Consequently, the present method exhibits a considerable 
advantage in time savings over finite element analyses in a study of this nature. 

[ f±45k fc and ( 45) ^ laminates with circular cutouts .- As an example of the capabilities of the 
current analytical method, the extreme case of a unidirectional (45)24 laminate exhibiting 
shear-extension coupling is analyzed and the results are compared with those for a more 
practical [(±45)6]s angle-ply laminate. The laminates contain circular cutouts and are 
subjected to a uniaxial compression load. The maximum normal stress, (o x ) m aximum. 
normalized by the applied compressive stress, g°, is plotted in figure 4 as a function of panel 
aspect ratio, L/W, for different cutout diameters. The ratio L/W varies from one to three, and 
D/W is equal to 0.1 , 0.3, and 0.6. The results for the [(±45>6]s laminate are represented by the 
solid line, the results for the (45)24 laminate are represented by the dashed line, and the 
symbols, shown at points where analytical values are calculated, represent the different 
cutout diameters. As the cutout size increases, the panel aspect ratio has a significant effect 
on the maximum normal stresses, which are much higher for low aspect ratio panels than for 
high aspect ratio panels. The maximum normal stresses increase for large cutout sizes, and 
maximum normal stresses in the (45)24 laminates are greater than maximum normal stresses 
in the [(±45)6]s laminates. 

The analytical results for the [(±45)6]s and (45)24 laminates subjected to an applied 
shear stress are shown in figure 5. The maximum shear stress, (o xy ) ma ximum. normalized by 
the applied shear stress, o xy , is plotted as a function of panel aspect ratio, L/W, for D/W equal 
to 0.1 , 0.3, and 0.6. Panel aspect ratio does not significantly influence the magnitude of the 
maximum shear stresses in a panel with a small cutout but becomes increasingly important 
as the cutout diameter grows larger. As with the compression-loaded laminates, the (45)24 
laminate has higher maximum shear stresses than the [(±45)e]s laminate under shear 
loading for the range of panel aspect ratios investigated in the present study. 

Analytical results are also calculated for the [(±45)6]s and (45)24 laminates subjected to 
combined compression and shear. The maximum normal stress, (o x )maximum. normalized by 
the applied compressive stress, o°, is plotted in figure 6; and the maximum shear stress, 
(o xy )maximum. normalized by the applied shear stress, o xy , is plotted in figure 7. The maximum 
shear and normal stresses are plotted as functions of panel aspect ratio. An examination of 
the results for the laminates loaded in shear only and the laminates subjected to combined 
shear and compression loading reveals that the addition of the compression load to the 
shear-loaded panels does not greatly impact the trends and magnitudes of the maximum 
shear stresses. However, for combined loading conditions, the maximum normal stresses 
exhibit behavior radically different from the maximum normal stresses in the laminates 
loaded in compression only (see fig. 4). Under combined shear and compression loading, 
the maximum normal stresses in the (45)24 laminates with D/W equal to 0.1 are lower than 
the maximum normal stresses in the [(±45)6]s laminates with the same D/W. For large cutouts 
and low panel aspect ratios, the maximum normal stresses in the (45)24 laminates are larger 
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than the maximum normal stresses in the [(±45)6]s laminates. However, for large cutouts and 
high aspect ratios, the maximum normal stresses in the [(±45)6]s laminates are higher than 
the maximum normal stresses in the (45)24 laminates. 

f (±45/0/90h k . I (±451Q2 ) 3 ]s , and f ( ±45/90 ? U k laminates with circular cutouts .- In order to study 
the effects of panel aspect ratio and cutout size on the stress field of some commonly used 
laminates, [(±45/0/90)3]s, [(±45/02)3]s, and [(±45/902)3]s laminates subjected to compression, 
shear, and combined loads are analyzed. The maximum normal stress, (a x ) ma ximum. 
normalized by the applied compressive stress, o x , is plotted in figure 8 as a function of panel 
aspect ratio for these three laminates subjected to a compression load for D/W equal to 0.1 , 
0.3, and 0.6. The 0° plies have the highest extensional modulus and, therefore, the 
[(±45/02)s]s laminate, which contains the most 0° plies of the laminates studied, has the 
highest maximum normal stresses for all cutout diameters. The [(±45/902)3]s laminate, whici 
does not contain any 0° plies, has the lowest maximum normal stresses for all values of DM'. 
The panel aspect ratio, L/W, has a minimal effect on the results for laminates with D/W equa 
to 0.1 . However, as the cutout diameter increases, the laminates with larger cutout diameters 
have very high maximum normal stresses at low panel aspect ratios, and these maximum 
normal stress values tend to coalesce as the panel aspect ratio increases. 

The maximum shear stress, (o X y) ma ximum. normalized by the applied shear stress, a° y , 
is plotted in figure 9 as a function of panel aspect ratio for the [(±45/0/90)3]s, [(±45/02)3ls. and 
[(±45/902)3 ]s laminates subjected to a shear load. For cutout-diameter-to-panel-width ratios, 
D/W, of 0.1 and 0.3, the [(±45/902)3]s laminate has the highest maximum shear stresses and 
the [(±45/0/90)3]s laminate has the lowest maximum shear stresses for the range of L/W 
considered. As the panel aspect ratio increases, this trend remains unchanged for D/W 
equal to 0.1 and 0.3. However, for D/W equal to 0.6, the ordering of the [(±45/0/90)3]s, 
[(±45/0 2 ) 3 ]s, and [(±45/902)3]s laminates in terms of the highest and lowest maximum shear 
stresses changes as the panel aspect ratio increases. For low values of L/W, the same trend 
that is observed for laminates with D/W equal to 0.1 and 0.3 occurs for D/W equal to 0.6. As 
the panel aspect ratio begins to increase, the [(±45/902)3]s laminate continues to have the 
highest maximum shear stresses, but the [(±45/02)3]s laminate has the lowest maximum 
shear stresses. For a panel aspect ratio of three, the [(±45/02)3]s laminate has the highest 
maximum shear stress, and the [(±45/0/90)3]s laminate has the lowest maximum shear stress. 

Results are also obtained for the [(±45/0/90)3]s, [(±45/02)3]s. and [(±45/902)3]s 
laminates subjected to combined shear and compression loads. The maximum normal 
stress, (a x ) ma ximum. normalized by the applied compressive stress, a°, is plotted as a function 
of panel aspect ratio in figure 10. The addition of the shear load to the compression load 
increases the maximum normal stresses. For all values of D/W, the [(±45/02)3]s laminate still 
has the highest maximum normal stresses, and the [(±45/902)3]s laminate still has the lowest 
maximum normal stresses. However, the change in the magnitude of the maximum normal 
stresses as the panel aspect ratio increases is not as large for the combined load case as for 
the compression load case. Maximum shear stresses, (a xy ) max j mum , normalized by the 
applied shear stress, a xy , are plotted as a function of panel aspect ratio in figure 1 1 . The 
addition of the compression load to the shear load increases the maximum shear stresses in 
the laminates, but the general trends remain the same with one exception. For panel aspect 
ratios near 1.0, the [(±45/02)3]s laminate has higher maximum shear stresses than the 
[(±45/902)3]s laminate for D/W equal to 0.3 and 0.6. As expected from symmetry 
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considerations, the [(±45/902)3]s and [(±45/02)3]s laminates have identical maximum shear 
stress values for a panel aspect ratio of 1 .0 when they are subjected to shear loading only 
(see fig. 9). 

Annle-olv laminates with circular cutouts .- Angle-ply laminates with circular cutouts are 
analyzed in order to examine the effects of fiber angle orientation on the maximum stresses 
in panels subjected to shear, compression, and combined loads. In figure 12, maximum 
normal stresses, (o x ) ma ximum. normalized by the applied compressive stress, <j°, are plotted as 
a function of fiber angle, 9, for square, compression-loaded, angle-ply laminates with circular 
cutouts. The solid lines represent analytical results for D/W equal to 0.1 , 0.3, and 0.6. The 
fiber angle, 0, for the [(±0) 6 ]s angle-ply laminate ranges from 0° to 90°. The maximum normal 
stresses are largest for a fiber angle of 0°, which coincides with the direction of the applied 
stress. The maximum normal stresses decrease steadily as the fiber angle is increased. This 
reduction in maximum normal stress values is attributed to the decrease in the extensional 
modulus of the angle-ply laminate as the fiber angle increases. 

The maximum shear stresses, (Oxy)maximumi normalized by the applied shear stress, 
c° r are plotted in figure 13 as a function of fiber angle, 0, for square, shear-loaded, angle-ply 
laminates with circular cutouts. The ratio D/W is equal to 0.1 , 0.3, and 0.6, and 0 ranges from 
0° to 90°. In some cases, the exact maximum stress is probably not found because the 
laminate stresses are calculated at a finite number of points on the panel. This discretization 
problem can cause some waviness in the maximum stress curves. The largest maximum 
shear stresses occur for a fiber angle of about 45° and decrease to minimums at fiber angles 
of 0° and 90°. This behavior appears to be consistent with the decrease in the laminate 
shear modulus as the fiber angle changes from 45°. 

Square angle-ply laminates with circular cutouts with D/W equal to 0.1 , 0.3, and 0.6 
are also analyzed for combined shear and compression loads. The maximum normal stress, 
(ox)maximum, normalized by the applied compressive stress, a°, and the maximum shear 
stress, (oxy)maximum, normalized by the applied shear stress, o° y , are plotted as functions of 
fiber angle in figures 14 and 15, respectively. The maximum normal and shear stresses 
show significant increases in magnitude over the maximum normal and shear stresses for the 
single load cases. The highest maximum normal stresses occur at a fiber angle of 0° as they 
did in the laminates subjected to a compression load only. The highest maximum shear 
stress for laminates with cutout-diameter-to-panel-width ratios, D/W, equal to 0.1 occurs at a 
45° fiber angle and shifts toward a fiber angle of 40° as D/W increases to 0.6. 

Angle-plv laminates with elliptical cutouts .- Angle-ply laminates with elliptical cutouts are 
examined to assess the effects of cutout shape, size, and orientation on the maximum 
stresses. For this analysis, the ratio of the major axis of the ellipse, D, to the panel width, W, 
is given values of 0.1 , 0.3, and 0.6. The ratio of the minor axis, A, to the major axis of the 
ellipse, D, is equal to 0.5, 0.75, and 1 .0. The major axis of the elliptical cutout is inclined at an 
angle 'P to the y-axis, and *¥ ranges from 0° to 90°. The maximum normal stress, (c x ) max j mum , 
normalized by the applied compressive stress, a°, is plotted as a function of *¥ for square, 
compression-loaded [(±30) 6 ]s. [(±45) 6 ]s. and [(±60) e ]s laminates in figures 16, 17, and 18, 
respectively. Similar to the results for the angle-ply laminates with circular cutouts, the 
maximum normal stresses are highest in the [(±30)6]s laminate and decrease as the fiber 
angle increases. Laminates having elliptical cutouts inclined at small angles to the y-axis 
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have the largest maximum normal stresses, and, in most cases, these maximum normal 
stresses increase as A/D decreases. As *F is increased, the maximum normal stresses 
decrease in all of the laminates with elliptical cutouts. At large values of *F the maximum 
normal stresses decrease as A/D decreases. 

Maximum shear stresses, (a X y) ma ximum. for the [(±30) 6 ] s , [(±45) 6 ]s, and [(±60) 6 ]s 
laminates subjected to shear loads are normalized by the applied shear stress, a x , and 
plotted as a function of the elliptical cutout orientation angle, 'F, in figures 19, 20, and 21 , 
respectively. The normalized maximum shear stress curves behave in a somewhat erratic 
manner. As mentioned previously, the laminate stresses are surveyed at a finite number of 
points. Therefore, the exact maximum stress is probably not always found, and some 
waviness in the maximum stress curves results. This waviness is magnified by the sensitivity 
of the value and location of the maximum shear stress to changes in V P, especially at low 
values of A/D. An example of this phenomenon is shown for the [(±45)6]s laminate in figures 
22 through 25. The shear stress distribution and the maximum shear stress location for 
laminates with the major axis of the elliptical cutout inclined at 0°, 25°, and 45° are shown in 
figures 22, 23, and 24, respectively. In contrast, for circular cutouts the maximum shear stress 
location remains fixed. For example, the shear stress distribution for the [(±45) 6 ]s laminate 
with a circular cutout is shown in figure 25. 

Although the maximum shear stress curves do not vary smoothly, the general trends 
are discernable. In the case of the [(±30)6]s laminate, the maximum shear stresses occur 
when *F is about 30° for the smaller cutouts and about 45° for the largest cutout size. In the 
[(±45)6]s laminate, the maximum shear stresses occur when 'F is equal to 45° for all cutout 
sizes. Under shear loading conditions, the maximum shear stresses in the [(±60)e]s laminate 
occur when *F is about 60° for small cutouts and about 45° for the largest cutout size. Of the 
three angle-ply laminates examined, the [(±45)6]s laminate has the overall highest maximum 
shear stresses. In all of the angle-ply laminates with cutout-major-axis-to-panel-width ratios, 
D/W, equal to 0.1 and 0.3, the maximum shear stresses increase as A/D decreases. For the 
laminates with D/W equal to 0.6, the maximum shear stress decreases as A/D decreases for 
values of *P near 0° and 90°. For some range of values of *F between 5° and 85°, depending 
on the particular laminate, the maximum shear stresses are highest for A/D equal to 0.5, and 
lowest for A/D equal to 0.75, with the maximum shear stresses for A/D equal to 1 .0 falling 
between the two. 

The maximum normal stresses, (a x ) maX imum, normalized by the applied compressive 
stress, o°, are plotted as a function of ¥ for the [(±30) 6 ]s, [(±45) 6 ]s, and [(±60) 6 ]s laminates 
subjected to a combined shear and compression load in figures 26, 27, and 28, respectively. 
The addition of the shear load to the compression load causes the maximum normal stresses 
to increase to values approximately double the maximum normal stresses due to the 
compression load alone. In addition, the highest maximum normal stresses no longer occur 
near *F equal to 0°. Instead, they occur for values of 'F ranging from 20° to 40°. Similar to the 
compression-loaded laminates, the [(±30)6]s laminate has the highest maximum normal 
stresses; and, in most cases, the maximum normal stresses in the laminates under combined 
loading increase as A/D decreases. The maximum shear stresses, (a xy ) max imum, normalized 
by the applied shear stress, a° y , for the laminates subjected to combined loads are shown as 
a function of *F in figures 29, 30, and 31 . The addition of the compression load to the shear 
load increases the maximum shear stress values and changes the elliptical cutout rotation 
angle at which the highest maximum shear stresses occur. Similar to the shear-load-only 
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case, the [(±45)6]s laminate has the highest stresses, and a decrease in A/D results in an 
increase in the maximum shear stress. 


CONCLUDING REMARKS 


A method which combines Lekhnitskii's complex variable equations with boundary 
collocation and a Laurent series approximation is used to analyze the stress distributions in 
some finite anisotropic panels with centrally located cutouts. The maximum normal and 
shear stresses are found for panels with circular and elliptical cutouts and subjected to 
compression, shear, and combined loads. The effects of panel aspect ratio, anisotropy, 
cutout size, and cutout orientation on the magnitude and distribution of the stresses are 
examined. Angle-ply laminates, unidirectional off-axis laminates, as well as [(±45/0/90)3]s, 
[(±45/0 2 )3 ]s. and [(±45/902)3]s laminates are studied in order to demonstrate the capabilities 
of the analytical method. To assess the accuracy of the method, analytical results are 
compared with experimental data for orthotropic panels with circular cutouts. The panels are 
subject to either compression or shear loading. In both cases, the analytical results show 
good agreement with the experimental data. 

To demonstrate the flexibility of the analytical method with respect to material 
properties and panel geometry, several laminates subjected to compression, shear, and 
combined loads are analyzed. To study the effects of panel aspect ratio on laminate 
stresses, the limiting case of a (45) 2 4 laminate with a circular cutout is analyzed and 
compared to more common laminates constructed of 0°, 90°, and 45° plies. It is shown that 
the effect of panel aspect ratio on the magnitude of the stresses in panels with small cutouts 
is minimal, but as the cutout size increases, the effects become more significant. A study of 
square [(±0) 6 ] s angle-ply laminates with circular cutouts identifies relationships between 
shear and extensional moduli and the maximum normal and shear stresses. Angle-ply 
laminates with high extensional moduli, such as the ( 0) 2 4 laminate, have high maximum 
normal stresses when subjected to compression loading. As the fiber angle deviates from 0°, 
the maximum normal stresses decrease. Angle-ply laminates with high shear moduli, such 
as the [(±45)6]s laminate, have the highest maximum shear stresses when subjected to shear 
loading, and as the fiber angle deviates from 45°, the maximum shear stresses decrease. 

A study of square angle-ply laminates with elliptical cutouts whose major axes are 
inclined to the panel coordinate system is conducted to demonstrate the ability of the current 
analytical method to accommodate different cutout geometries. The results indicate that the 
values of the maximum stresses depend on the relationships between several geometric 
factors and material properties. For example, the ellipse minor-axis-to-major-axis ratio has 
an effect on the maximum stresses that depends on the cutout size, the fiber angle of the 
angle-ply laminate, and the angle of inclination of the elliptical cutout. As the inclination of 
the elliptical cutout is varied, the stress distribution changes dramatically, especially for shear 
loading where the location of the maximum shear stress also changes. 

The wide range of results obtained in the present study exemplifies the ability of the 
present method to adapt to many different panel sizes and cutout geometries. The method 
allows for biaxial, shear, and combined loads. While only traction loading conditions are 
considered in the present paper, the method is also capable of analyzing displacement 
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loading conditions. The flexibility of this method characterizes its main advantage over finite 
elements in the areas of structural design and optimization. In addition, the FORTRAN 
program that implements the present method requires minimal computer time and can run 
interactively. Because the present method requires only two lines of input to describe the 
model, user modeling time is significantly less than with finite element analyses. 
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Fig. 1. 


Fig. 2. 


Table 1 . Lamina properties of Hercules, Inc. AS4/3502 graphite-epoxy 
material 


Longitudinal Young's Modulus, psi 
Transverse Young's Modulus, psi 
Shear Modulus, psi 
Major Poisson's ratio 


18.5 x 10 6 
1.6 x 10 6 
0.832 x 10 6 
0.35 



Analytical model of a shear- and compression-loaded rectangular plate with a centrally 
located elliptical cutout inclined to the y-axis. 

x 



1 1 1 1 i i 

0 0.2 0.4 0.6 0.8 1.0 

W/2 

Stress distribution near a circular cutout in a rectangular (±45/02/±45/02/±45/0/90)2s 
panel under applied uniaxial compression (experimental data from reference 6). 
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Analysis 



Fig. 4. Maximum normal stresses in rectangular uniaxial-compression-loaded [(±45) 6 ]s and Fig. 6. Maximum normal stresses in rectangular compression- and shear-loaded f(±45)els and 
(45)24 laminates with a circular cutout. (45) 24 laminates with a circular cutout. 






D/W [[(-45/0/90 ) 3 ] s [(-45/0 2 ) 3 ] s [(-45/90 2 ) 3 J 
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Fig. 12. Maximum normal stresses in square compression-loaded angle-ply laminates with a Fl 9 14 Maximum normal stresses in square compression- and shear-loaded angle-ply 
circular cutout. laminates with a circular cutout. 
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Fig. 24 Shear stress distribution for a square shear-loaded [(±45) 6 ]$ laminate with an Fig 26. Maximum normal stresses in a square compression- and shear-loaded [(±30)6]s 
elliptical cutout inclined at 45° to the y-axis. laminate with an elliptical cutout inclined to the y-axis. 
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Fig. 28. Maximum normal stresses in a square compression- and shear-loaded [(±60)e)s Fig. 30. Maximum shear stresses in a square compression- and shear-loaded [(±45)e]s 
laminate with an elliptical cutout inclined to the y-axis. laminate with an elliptical cutout inclined to the y-axis. 








X F, degrees 

Fig. 31 . Maximum shear stresses in a square compression- and shear-loaded [(±60)e]s 
laminate with an elliptical cutout inclined to the y-axis. 
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